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^ . Abstract 

We consider the weak decays of a B meson to final states that are mixtures of 
S-wave radially excited components. We consider non leptonic decays of the type 
B p'tt/B pn, B lu'tt/B — > ujtt and B (fy'ir/B (j)7r where p' , uj' and 0' are 
higher p, lo and (p resonances. We find such decays to have larger or similar branching 
ratios compared to decays where the final state p, uj and (p s-i's in the ground state. 
We also study the eff'ect of radial mixing in the vector and the pseudoscalar systems 
generated from hyperfine interaction and the annihilation term. We find the effects 
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of radial mixing to be small and generally negligible for all practical purposes in the 
vector system. However, in the r]—r}' system the effects of radial mixing are appreciable 
and seriously affect decay branching ratios for B — *■ •q{r}')K{K*). In particular we find 
that nonstandard rjij]') mixing can resolve the puzzles in S — > r]{r]')K decays. 
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1 Introduction 



Nonleptonic B decays play a very important role in the study of CP violation. It is expected 
that these studies will test the standard model(SM) picture of CP violation or provide hints 
for new physics. Most studies of two body nonleptonic B decays have concentrated on 
processes of the type B — > M1M2 where both Mi and M2 are mesons in the ground state 
configuration. Here we want to look at nonleptonic B decays to final states where one of the 
final state meson contain admixtures of radially excited components. We expect such decays 
to have larger or similar branching ratios compared to decays where the final state contains 
the same meson in the ground state. There is an easy explanation for such a statement. 
For simplicity let us consider a simple model in which B — > ttM and M is a simple flavor 
eigcnstatc with no flavor mixing beyond isospin; e.g. p, uj or (f) and we are interested in 
comparing the branching ratios for the final states containing the ground state configuration 
of meson M and its radial excitation. Some possible examples are: 

BR{B° ^ 7r-p+')/5i?(S° ^ 7r-p+) (1) 

BR(B- n-u') /BR{B- n-u) (2) 

BR(Bs TT V) /BR(B, ttV) (3) 

where p"*"', u' and 0' are radially excited states. 

We assume an extreme factorization approximation in which the b quark decays into a 
pion and a w-quark and we neglect the relative Fermi momentum of the initial b and the 
spectator q. The quark transition for the processes (1) and (2) is then 

b^7r-{p)u{-p) (4) 

where p denotes the final momentum of the tt" . For the process (3) the quark transition is 
essentially similar to the one above 

b ^ 7r'{p)s{-p) (5) 

where now p denotes the final momentum of the 7r°. 

Concentrating on the processes (1) and (2) the transition matrix for the full decay has 
the form 

{n-{p)M{-p)\T\B) = {M{-p)\F\u{-mO))-{^-{p)u{-p))\w\b) (6) 
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where T denotes the transition matrix for the hadronic decay which factors into a weak 
matrix element at the quark level denoted by W and a fragmentation matrix element denoted 
by F describing the transition of a quark with momentum —p and an antiquark with zero 
momentum to make a meson with momentum —p. 

It is immediately clear that if the final momentum p is large, the fragmentation matrix 
element will depend upon the high momentum tail of the meson wave function. This might 
tend to favor radial excitations over ground states, since the radial excitations are expected 
to have higher kinetic energies. We now note that the harmonic oscillator wave functions 
commonly used in hadron spectroscopy have a Gaussian tail for their high momentum com- 
ponents and this can suppress the fragmentation matrix element {M{—p) \ F p)g(O)) in 
comparison with wave functions from a different confining potential. Hence the branching 
ratio to a final state which is radially excited will be sensitive to the choice of the confining 
potential. 

So far we have assumed the physical states to be pure radial excitations. However ad- 
ditional interactions can mix the various radial excited components. For instance hyperfine 
interactions can mix radial excitations with the same flavor structure and so in general in the 
p, uj and (f) system the various physical states will be admixtures of radial excitations [|l], 0. 
Flavor mixing in the vector system is known to be small but is important in the pseudoscalar 
sector. Here the mixing in the 77(^7') system receives an additional significant contribution 
from the annihilation diagram that leads to flavor mixing of the strange and non strange 
parts of the 77(77') wavefunction. The mixing in the pseudoscalar sector, therefore, is different 
from the ideal mixing found in the uo — cf) system. It is also possible for the annihilation term 
to mix states that are radial excitations allowing the the rjij]') wavefunction to contain radi- 
ally excited components . Such non standard 77(77') mixing can have important implication 
for the non leptonic decays B — > r]{ri')K{K*). 

In the transitions chosen in (1),(2) and (3) the radially excited meson must include 
the spectator quark and therefore must depend upon the fragmentation matrix element 
{M{—p)\F\u{—p)q{0)). For the case of B ^ rj{rj')K decays, factorization results in the 
kaon leaving the weak vertex with its full momentum and the remaining quark carries the 
full momentum of the final meson. There is therefore a form factor in which there is a large 
internal momentum transfer needed to hadronize this quark with the spectator antiquark. 
This might favor the radial excitation if it has a much higher mean internal momentum. 
However things are more complicated here as the 77(7/') can also be produced by an qq pair 
in a penguin diagram without containing the spectator quark. One possibility is when the s 
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quark in the QCD penguin combines with the s quark from the b ^ s transition to form the 
T]{ri'). Another possibihty is when a q and q pair (where q = u,d,s), appearing in the same 
current in the effective Hamiltonian, hadronizes to the Tjij]'). In the diagrammatic language 
this is often represented as a "gluon" sphtting into a qq pair which then hadronizes into a 
i]{ri'). This term is usually called OZI suppressed ^ as in most decays the contribution 
of this term is indeed suppressed with respect to other terms in the decay amplitude. This 
may not be the case in the B —* r]{r]')K decays where the OZI term can be of comparable 
size as other terms in the decay amplitude and in particular we show that B — > t]'K* can 
have significant contribution from the OZI suppressed term. 

This paper is organized in the following manner: In the next section we study the mass 
mixing in the vector meson sector, involving the p, u and 0, and the pseudoscalar sector 
involving the ri{ri'). In section 3 we present a general treatment of nonleptonic B decays 
using the effective Hamiltonian and the factorization assumption. We then show how this 
approach is related to the diagrammatic approach of studying nonleptonic B decays. In 
section 4 we study nonleptonic decays of B to final states involving higher resonant p, u 
and states. This is followed by section 5 where we make predictions for B ri{ri')K and 
B r]{ri')K* and comment on the relevance of the OZI suppressed term in the calculation 
of these decays. Finally we present our conclusions. 



2 Mass mixing in the vector meson and pseudoscalar 
sector 

We start with the mixing for the p system. To obtain the eigenstates and eigenvalues we 
diagonalize the mass matrix which has the form 

< %%,n\M\qaqb,n > = ^aa'hw^nn'inia + mfe + En) + 6aa'hb' Sa " S6?/'„(0)V^„/ (0)(7) 

maTrib 

where and nia^h are the quark spin operators and masses. Here n = 0, 1, 2 and the basis 
states are chosen as >= \uu + dd > / \/2 and S* = |ss >. In the above equation En is 
the binding energy of the n^^ radially excited state and B is the strength of the hyperfine 
interaction. Note we are only presenting results for the neutral mesons. A similar treatment 
also can be used for the charged mesons. We will use a very simple model for confinement in 
our calculations as we do not intend to present a detailed study of light meson spectroscopy. 
Our aim, as already stated in the previous section, is to study non leptonic decays of the B 
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meson to radially excited light meson states as well as to study the effects of radial mixings 
in the non leptonic decays B — * T]{ri')K{K*). We believe the conclusions reached on the 
basis of our calculations are likely to hold true in a more detailed model of confinement. 

To begin with, we use the same harmonic confining potential as well as the other param- 
eters used in Ref0 to obtain the eigenstates and eigenvalues for the mass matrix in Eqn. ^. 
The various parameters used in the calculation are m„ = = 0.350 GeV, = 0.503 GeV, 
the angular frequency, u =0.365 GeV and b = B /m^ =0.09. 

We obtain for the eigenvalues and eigenstates in the p system 

|p(0.768) > = 0.990|Ar >o +0.124|iV >i -0.066|Ar >2 
|p(1.545) > = 0.108|A^ >o -0.973|A^ >i +0.204|Ar >2 
|p(2.370) > = 0.089|A^ >o -0.195|iV >i +0.977|A^ >2 (8) 

To see how this result changes with a different confining potential we use a power law 
potential V{r) = Ar" P| . We will use a linear and a quartic confining potential and 
compare the spectrum with that obtained with a harmonic oscillator potential. To fix the 
coefficient A we require that the energy eigenvalues of the Schrodinger equation are similar 
in the least square sense with the energy eigenvalues used in Ref0]. So for example, for the 
linear potential, we demand that 

F = y^{En{harmonic) — En{linear)Y 

n 

is a minimum. This fixes the constant A in V{r) = Ar and we obtain 

|p(0.775) > = 0.992|Ar >o +0.112|iV >i -0.053|Ar >2 
|p(1.515) > = 0.104|iV >o -0.986|iV >i +0.130|iV >2 
|p(2.260) > = 0.066|A^ >o -0.122|iV >i +0.99|A^ >2 (9) 

We follow the same procedure for the quartic potential and obtain 

|p(0.759) > = 0.988|A^ >o +0.129|iV >i -0.077|A^ >2 
|p(1.567) > = 0.103|iV >o -0.955|iV >i +0.278|A^ >2 
|p(2.550) > = 0.11|A^ >o -0.267|iV >i +0.957|A^ >2 (10) 

We observe that the mass eigenstates and eigenvalues are not very sensitive to the con- 
fining potential and the radial mixing effects are small. 
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We next turn to mixing in the t] — rj' system. In the traditional picture the t] and t]' 
mesons are mixtures of singlet and octet states rji and rjs of SU{3). 

cos 6 — sin 9 
sin 9 cos 9 

1 

uu + dd + ss 



V 




7]' 





^8 



UU + dd — 2ss 



1 

73 



(12) 
(13) 



where the mixing angle 9 lies between —10*^ and — 20°[|T( 



To obtain the eigenstates and eigenvalues in the r] — rj' system we diagonalize the mass 
matrix 

B 



< q'^ql n'\M\qaqb, n> = daa'^bb'^nn^rria + + E„) + daa'^i 



bb' 



rriamh 



+ ^ab^a'b' " 



A 



, , , , (14) 

rriamb 

This has a similar structure as the p system but now we have the additional annihilation 
contribution with strength A that causes flavor mixing. In our calculations we use the phase 
convention in Ref[@] where the wavefunctions at the origin in configuration space, which 
enter in the hyperfine and annihilation terms in the mass matrix, are positive (negative) for 
the even(odd) radial excitations. 

We try to fit the values of A and B to the measured masses. The mass matrix is 6 x 6 
matrix which we diagonalize to make predictions for 6 masses and mixings. However, for 
the sake of brevity we will only give the predictions for rj and 1]' masses and wavefunctions. 
Several solutions that give acceptable values of the masses can be obtained. We choose solu- 
tions for the linear, quadratic and quartic confining potentials that make similar predictions 
for the T]{ri') masses: 

For the linear potential we obtain with B = 0.065m^ and A = 0.045m^ 

|?7(0.544)) = 0.961|A^ >o -0.198|A^ >i +0.108|A^ >2 

- 0.150|S >o +0.050|^ >i -0.032|5 >2 

|V(0.924)) = 0.170|iV >o +0.0385|A^ >i -0.0157|A^ >2 

+ 0.974|S >o -0.126|S >i +0.0486|S >2 (15) 

For the harmonic potential we obtain with B = 0.065m^ and A = 0.065m^ 
|r7(0.547)) = 0.913|iV >o -0.252|A^ >i +0.154|A^ >2 
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- 0.249|^>o +0.107|5>i -0.076|S>2 
|r7'(0.931)) = 0.316|A^ >o +0.109|iV >i -0.049|Ar >2 
+ 0.92515 >o -0.1481^ >i +0.0881^ >2 



(16) 



Our results for the harmonic potential is similar to the results obtained in Ref|0] where a 
slightly different mass mixing matrix than the one used here has been used to obtain the 
1] — 1]' mixing. 

Finally, for the quartic potential we obtain with B = 0.065m^ and A = 0.1 Im^ 

|r/(0.547) > = 0.764|A^ >o -0.287|A^ >i +0.198|A^ >2 

- 0.441|5 >o +0.248|S >i -0.196|S >2 

|r/'(0.940) > = 0.623|A^ >o +0.350|iV >i -0.177|iV >2 

+ 0.658|5 >o -0.134|5 >i +0.087|5 >2 (17) 

It is clear that the eigenstates of ri{ri') system are sensitive to the confining potential and 
there can be substantial radial mixing which can then affect the predictions for the decays 
B ri{ri')K{K*) . We note that as we move from the linear to the quartic potential the 
i]{r]') mixing deviates more significantly from the ideal mixing case. This may be understood 
from the fact that the fit to annihilation term, A, is the largest for the quartic potential 
which leads to the largest deviation from the ideal mixing case. A standard t] — rj' mixing 
often used in the literature is given by ||, 0] 

1 
76 



\v')std 
\v)std 



1 

73 



uu + dd + 2ss 
uu + dd — ss 



We can then write the 77 — 77' states obtained with the various confining potential and keeping 
only the ground states, in terms of the the states defined above. For the linear potential we 
find 



For the harmonic potential 



W) = 0.89 |r/'),,, - 0.43 |r/),,. 
It/) = 0.87 |7/),,, + 0.43 It?'),,, 



W) = 0.94 |t^'),,,- 0.27 |t^),,. 
It/) = 0.89 It/),,, + 0.32 It/'),,, 



(19) 



(20) 
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and finally for the quartic potential one finds 

W) = 0.90 IV),,, + 0.13 |r^),,, 

\V) = 0.88 |r/),,, + 0.08 IV),,, (21) 

This shows that all three confining potentials give mixings for the rj — i]' that have substantial 
overlap with the standard mixing but the mixing with the quartic potential is closest to the 
standard mixing in the sense that here one has the smallest component of the l^]),,, i\v')std 
in IV) (|r/)). 

We now turn to the u — (f) system. As in the the t] — rj' system we diagonalize the mass 
matrix in Eqn. |14[ We use the same value for the hyperfine interaction as used for the p 
system. Again, for the sake of brevity, we only give the wavefunctions for the ground and 
the first excited states. 

For the linear potential we obtain with B = 0.09m^ and A = 0.005m^ 

|u;(0.782)) = 0.991|iV >o +0.123|A^ >i -0.058|iV >2 

- 0.014|5 >o +0.004|S >i -0.002|S >2 
10(1.05)) = 0.012|A^ >o +0.011|iV >i -0.004|A^ >2 

+ 0.997|5 >o +0.071|5 >i -0.034|S >2 

|tu(1.52)) = -0.113|A^ >o +0.982|iV >i +0.144|A^ >2 

- 0.006|5 >o -0.034|S >i +0.004|S >2 
10(1.66)) = 0.007|A^ >o -0.030|A^ >i -0.014|A^ >2 

+ 0.068|5 >o -0.994|5 >i -0.077|5 >2 (22) 

For the harmonic potential we obtain with B = 0.09m^ and A = O.OlSm^ 

|cu(0.783)) = 0.984|iV >o +0.154|iV >i -0.081|Ar >2 

- 0.033|5 >o +0.011|5 >i -0.007|S >2 
10(1.05)) = 0.026|iV >o +0.029|iV >i -0.011|iV >2 

+ 0.994|5 >o +0.089|5 >i -0.048|S >2 

|c<;(1.57)) = -0.126|A^ >o +0.948|A^ >i +0.256|iV >2 

- 0.008|5 >o -0.139|5 >i +0.010|5 >2 
10(1.68)) = 0.025|iV >o -0.118|A^ >i -0.07|iV >2 

+ 0.082|5 >o -0.976|5 >i -0.143|5 >2 (23) 
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For the quartic potential we obtain with B = 0.09m^ and A = 0.023m^ 

|cu(0.782)) = 0.980|iV >o +0.163|A^ >i -0.096|A^ >2 

- 0.049|5 >o +0.012|5 >i -0.009|5 >2 
\(j){l.05)) = 0.041|A^ >o +0.034|A^ >i -0.015|iV >2 

+ 0.991|5 >o +0.100|5 >i -0.060|S >2 

|cu(1.58)) = -0.122|iV >o +0.932|A^ >i +0.322|iV >2 

- 0.010|5 >o -0.113|5 >i +0.006|5 >2 
10(1.7)) = 0.022|A^ >o -0.089|A^ >i -0.067|A^ >2 

+ 0.086|5 >o -0.968|5 >i -0.207|5 >2 (24) 

As in the p system we find the mixing to be insensitive to the confining potential and we 
also find the effects of radial mixing to be small. We also find, as expected, a smaller value 
for the annihilation term in the fits to the masses as compared to the pseudoscalar system. 



3 Effective Hamiltonian, Factorization and the Dia- 
grammatic approach 

In the Standard Model (SM) the amplitudes for hadronic B decays are generated by the 
following effective Hamiltonian p|: 

^ 10 

H!ff = -^[Vf,v;^ic,oif + c,oi^) - + + v,,v:/;)Oi] + H.c.{2h) 

i=3 

where the superscript m, c, t indicates the internal quark, / can be m or c quark and q can 
be either a c? or a s quark depending on whether the decay is a AS* = or AS = — 1 process. 
The operators Of are defined as 

Olf = q^^.Lfpf^YLh^ Olf = qj,Lfh^Lb 

01,5 = QliiLbq-ff,L{R)q' O^q = qa,-f^,Lhpqp-f^L{R)q'^ (26) 

3 3 
^7,9 = ^qif,Lbeg:qYRiL)q 0|io = l^Qali,Lhpeq>qp-i^R{L)q^ 

where R{L) = 1 ± 75, and q' is summed over u, d, and s. Oi and O2 are the tree level 
and QCD corrected operators. O^-q are the strong gluon induced penguin operators, and 
operators Oj-io are due to 7 and Z exchange (electroweak penguins), and "box" diagrams at 
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loop level. The Wilson coefficients c{ are defined at the scale /i ~ nib and have been evaluated 
to next-to- leading order in QCD. The c* are the regularization scheme independent values 
obtained in Ref. 0. We give the non-zero c{ below for nit = 176 GeV, as{niz) = 0.117, and 
fi = nif, = 5 GeV, 



ci = -0.307 , C2 = 1.147 , 4 = 0.017 , 4 = -0.037 , 4 = 0.010 , 4 = -0.045 

4 = -1.24 X 10"^ , 4 = 3.77 X 10~^ , 4 = -0.010 , 4o = 2-06 x 10"^ , 

U,C _U,C /7\r 1~>U.C / AT ^U,C T>U.C ,u,c 

3,5 — 



cl%lN, = PriNc. c^,% = P:'\ 4,1^0 = (27) 



where is the number of colors. The leading contributions to P]g are given by: PI = 
(§)c2(f +G'K,/i,g2))andP: = (^^)(iV,ci+C2)(f /x, g^)). The function G(m, /i, g^) 
is given by 

G'm,/i,g2 =4/ x 1 - x In \ ^ dx 28 

All the above coefficients are obtained up to one loop order in electroweak interactions. The 
momentum q is the momentum carried by the virtual gluon in the penguin diagram. When 
> 4m^, G{ni, fi, q^) develops an imaginary part. In our calculation, we use, for the current 
quark masses at the scale /i ~ m;,, m„ = 5 MeV, = 7 MeV, nig = 100 MeV, nic = 1.35 
GeV 0, [111 and use = m^/2. 

The structure of the effective Hamiltonian allows us to write the amplitude for B 

Ml Ms as 

Gt 

Now we can write the tree amplitude T as 

T = {Mi,M2\ciOif + C202f\B) 

T = C2{M^,M2\02f-{l/Nr)Oif\B) (30) 

where from Eqn. ^ A^i = — Cs/ci = 1.147/0.307 = 3.73. In the factorization assumption 
there are two contributions to the tree matrix element, T. To be specific consider the decay 
B- p^{uj)'n-. In this case there can be two contributions to T, given in the factorization 
assumption by 

T = T1+T2 

Ti = c2<Tx-\d^^Lh\B- >< p\u)\u^^Lu\Q>[^-^] 

T2 = C2<p'iuj)\uYLb\B' ><n~\d^^Lu\0>[l--^] (31) 
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In the diagrammatic language the first term, T2, corresponds to a 6 quark transition to a m 
quark and a W~ which turns into a vr^ . The u quark then combines with the spectator quark 
to form the p^{uj) particle. In the term Ti, the u quark from the h ^ u transition combines 
with the u quark from the W~ to form the p^{uj) particle while the d quark from the W~ 
combines with the spectator quark to form the tt~ . This is the color suppressed diagram and 
from the expression above we see that there is an additional suppression coming from the 
Wilson's coefficients and so, effectively Ti is suppressed by a factor of (1/A^c — V^i) ~ 1/15 
relative to T2. 

We now turn to the penguin contribution, and for simplicity, we just concentrate on the 
the t penguin, Pt. We can write 

Pt = {MuM2\J240i\B) (32) 

i 

Again from the values of the Wilson's coefficients given in Eqn. ^ we can write 

Pt ^ {M,,M2\clOs + clO,\B) + {M,,M2\cl05 + 40e\B) 

- ^4(Mi,M2|09|5) 

Pt ^ cUMi,M2\0^- {l/N2)0s\B) + cUMi,M2\0e- {l/N,)0,\B) 

- ^4(Mi,M2|09|5) (33) 

where from Eqn. |3 A^'s = -c*/4 = 0.037/0.017 = 2.2 and A^s = -4/4 = 0.045/0.010 = 4.5. 

In the diagrammatic approach the q'q' quarks appearing in the operator O3 — Oq appears 
from a "gluon" splitting into a q'q' pair while in the case the operators O7 — Oiq it is a "7" , 
"Z" splitting into a q'q' pair. 

Concentrating on only the term proportional to c\, we can write in the factorization 
assumption. 



Pt = Ptl + Pt2 

Pti = 4(1 ^ 



< p°{u)\u-f>'Lb\B- >< n-\djf,Lu\0 > 



+ < 'K-\d^^'Lb\B- >< p°(cj)|J7^Ld|0 > 
Pt2 = 4 < 7r-\dYLb\B- >< p'iLu)\u^^Lu + d^^Ld\0 > [-1 - -1] (34) 

iVc iV2 

We see that the second term, Pt2 has a suppression factor of {1/Nc — I/N2) ~ 1/8. This 
term is called OZI suppressed and in the diagrammatic language this is shown as a "gluon" 
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splitting into a quark-antiquark pair which then hadronizes to a hadron. Of course for a real 
gluon this process is forbidden as the color octet gluon cannot form a color singlet hadron. 
Note that there can be other OZI violating diagrams that have been considered to explain 
the large branching ratios in the decay B —>■ rj'K and the semi- inclusive decay B rj'Xg. 
In these diagrams the enhanced branching ratios are due to the anomaly, gluon couplings to 
the flavour singlet component of the rj' or the intrinsic charm content of the 1]' |jT2|, |13|. We 
will not consider such diagrams in our analysis. 

The terms represented by Pa, in the diagrammatic approach, has a "gluon" splitting into 
a quark- antiquark pair but now the antiquark combines with the d quark, coming from the 
b d transition, to form a meson while the other quark combines with the spectator quark 
to form the second meson in the final state. The two terms represented in Pfi represent the 
cases where the quark-antiquark pair from the "gluon" is a uu and a dd pair. 

One can do the same exercise with the term proportional to Cg in Eqn. ^ and in this 
case the OZI violating term is suppressed by {1/Nc — I/N3) ~ 1/9. Note the term Pt2 from 
the electroweak penguin term Cg does not have any suppression. This is expected as a qq 
pair from a 7 or Z boson is in a color singlet state and therefore can form a hadron. 

For the case of B ri{ri')K the term Pti and Pt2 are similar to the one above. 

Pt = Pn + Pt2 

Pti = c\{l-j;^j[<v{v)\uYLb\B-><K-\s-f,Lu\0> 

+ < K~\sYLb\B- >< r]{r]')\s-^f,Ls\0 > 
Pt2 = c\< K-\s^^'Lb\B- >< r]{r]')\u-f^Lu + dj^Ld + s-f^Ls\0 >[-*-- (35) 

Note we again find the OZI term to be suppressed by a factor ~ 1/8. However if terms in 
Pfi interfere destructively then the OZI terms may become important. Note the suppression 
in the OZI term can also be diluted if the contributions from the u,d and s term interfere 
constructively in Pt2- 

In the decay B^ p^{uj)'n-^ the OZI suppressed term does not play an important role, as 
these decays are not penguin dominated because the CKM factors in the tree and penguin 
terms are of the same order and the penguins are loop suppressed. This fact is also supported 
by recent experimental measurement of the branching ratios BR[B^ —>■ uj-k^] ~ 11.3 x 10^^ 



and BR[B —>■ p^n ] ~ 10.4 x 10 ^ |14|. Note that in Eqn. |3^, from the flavor structure of 



the and u wavefunction, it is obvious that the two terms in in Pn interfere destructively 
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for the but constructively for the uj. Furthermore the OZI suppressed term, Pt2 for the 
p° vanishes, neglecting the electroweak contribution, but not for the uo. This means that the 
penguin term for pvr^ is smaller than in uj-k^ . Hence if the penguin terms 

were dominant then there would be a significant difference between the branching ratios 
BR[B~ —>■ uJTf^] and BR[B^ P^r^]. Hence the small measured difference in the branching 
ratios for B^ —>■ p^Ti^ and B^ uj-k^ implies relatively small penguin effects in these decays. 

However decays of the type B — rjK are dominated by penguin terms because of large 
CKM factors in the penguin terms compared to the tree term. Here, the OZI suppressed 
terms may have significant effects on the predictions of these decays. Note that we expect 
the OZI suppressed terms to be more important in i? — > KP than in i? ^ KV decays where 
P is a pseudoscalar and \^ is a vector state. This follows from the fact that in J/ip and 
T decays we know that the OZI-forbidden process requires three gluons for coupling to a 
vector meson and two gluons for coupling to a pseudoscalar. Thus one would expect that 
the contribution of the OZI suppressed term should be much smaller in the B — > Kp^iuo) 
and B Kip decays than m B ^ Krj and B Krj' decays |]TB|, One of the authors 



of this work has shown that one can make definite predictions about the branching ratios 
B TjK/B Tj'K and B r]'K*/B rjK* ^ if one assumes that the OZI terms 
are forbidden. We will first derive these predictions in the language of effective Hamiltonian 
and using the factorization assumption. We then study how the predictions change if we use 
non standard mixing in the rj — r]' sector and if we include the OZI terms. 

4 B ^ p^uj^cf) transitions. 

In this section we study decays of the type B —>■ VP where V = p,uj,(j). As we found in 
section 2 the wavefunction of V has the general form, 

\V) = (36) 

i 

We can then write, 

Amp{B^VM) = Y.^^{M^N^\Heff\B) (37) 

i 

We now consider the ratios of the following decays 

Rp+ = BR{B° ^n-p+')/BR{B^ ^n-p+) (38) 
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RpO = BR{B- TT-p"')/BR{B- -n- p^^ 
R^ = BR{B- 7r-Lo')/BR{B- tt-lo) 
R^ = BR{Bs 'K°(j)')/BR{Bs iT^(j)) 



(39) 
(40) 
(41) 



As discussed in the previous section we can neglect the penguin contribution to these 
decays. Note that for the decay Bg TT^(f) there is no contribution from the QCD penguin 
and the dominant electroweak penguin term has the same structure as the tree amphtude 
and hence the ratio R^ remains essentially the same even in the presence of penguin terms. 
We then obtain 

Rn 









cJ7^(1-75)m|0) 






50) (tt- 


J7^(1-75)m|0) 



A' 



A, 



2 p3 



(42) 



where P is the magnitude of the three momentum of the final states and the form factor A^ 
is defined through 



{Vf\A^\P,) = {M, + Mf)A, 
+ 2MfAo^-^q^, 



e.q 



-A<, 



e.q 



f 



M? 



(43) 



RpO 







1^7^(1 


-75)b 


B^) 


> (vr 1 ^7^(1 


-75)m|0) 






U7'^(l 






(7r-|(i7^(l 


-75)w|0) 






2 p3 
r-po, 












Ao 













(44) 



In the above we have ignored the term in the amplitude for B p{p')tt which is given 
by Ti in Eqn. ^ and hence is suppressed by ~ 1/15 relative to the dominant term T2. 

Within the same approximation we can then write 



Ru 





M7^(l 




50) 


(tt M7m(1 


- 75)^ |o) 


{CO 


'U7'^(l 


-75)^ 




(7r-|J7^(l 


- 75)w |0) 














Ao 


p3 

UJ 











(45) 
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Finally, 



Rtt 







^^7m(1 -75)m|0) 


(01 57^^(1 -75)& 


so) (ttO 


^^7At(l-75)M|0) 



p3 
-'0/ 



(46) 



To calculate the ratios we need the form factor Ao. We use the model for form factors 



in Ref|T9[ that incorporates some relativistic features and is relatively simple to use. The 
calculation of the form factors in this model is not always in agreement with experimental 
results, however, we expect the model to make reliable predictions for the ratio of form 
factors. We believe this feature to be true for most other models for form factors. Hence we 
shall mostly calculate quantities that can be expressed as ratios of form factors. The model 
in |]19| assumes, like many other models, the weak binding limit which sets the meson mass 
to be the sum of the masses of the constituent quarks making up the meson. Hence the effect 
of binding are included in the quark masses. This is a good approximation for the lowest 
resonances of the p, uj and the (f) systems. In this model the form factor Aq is given by 



Ao 

A 
Jl 

J2 



A[Ji- 



-J, 



AMiMfM^ 



J ■' iia^ nif' 



where 



M. 



± 



Mf ± Mi 
2ms0 = 2 

M - 



Ml 
rriimf 



(47) 



rrii + ruf 

and (pf and (pi represent the momentum space wave functions while (3 is the velocity of the 
mesons in the equal velocity frame and rriij are the non spectator quark masses of the initial 
and final meson. 
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We use the momentum wavefunction 0/ obtained from spectroscopy in section. (2) while 
for 0j we use the wave function 

0i = 0B = iVse-^''/^'^ (48) 

where is the Fermi momentum of the B meson. In our calculations we will take Pf = 
SOOMeV. Note that in the analysis presented in the introduction we have neglected the Fermi 
momentum of the b quark, but since pp/^b is small the general conclusions reached in the 
introduction still continue to be valid. 

For transitions to higher resonant states, we use the same quark masses as those used 
in the transition of the B meson to the lowest resonant state. This is reasonable, as the 
spectator quark still comes from the B meson and therefore has the same value for its mass 
irrespective of whether the final state is in the lowest or the first excited state. The values 
for the masses of the non spectator masses are taken to be essentially the same as those used 
in section. (2) for spectroscopy. However for the calculation of the velocity (3 and hence a 
defined in Eqn. ^ we use the physical mass of the higher resonant state. 

Table 1: Ratios of branching ratios for different confining potentials 



Ratio 


Linear 


Quadratic 


Quart ic 


Rp+ 


2.3 


2.0 


1.9 


RpO 


2.3 


2.0 


1.9 


Ruj 


3.5 


2.5 


1.7 


R(j, 


6.7 


6.2 


5.2 



In Table. |I] we give our predictions for the various ratios defined above. We find that the 
transitions to higher excited states can be comparable or enhanced relative to the transitions 
to the ground state. From Table. we see that the ratios of branching ratios are slightly 
sensitive to the confining potential and the ratios of branching ratios increase as we go from 
the quartic to the linear potential. This is because the wavefunction for the linear potential 
has a longer tail and hence more high momentum components than the wavefunction for 
the quadratic and the quartic potentials. The wavefunction for the quadratic potential, has 
in turn, a longer tail and hence more high momentum components than the wavefunction 
for the quartic potential. As mentioned in section 1 the form factor m B ^ M transition, 
where M is a light meson, is sensitive to the high momentum tail of the meson wavefunction. 
Hence we would expect the hierarchy (Ai^o)«inear- > {Aifl) quadratic > (^i,o)gMariic where Ai^q 
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are the form factors for the transition of B to the first radially excited and the ground state 
of the meson M. We see from from Table. |I| that this hierarchy is maintained for the ratios 
of form factors and so we have {Ai/ Ao)iinear > {Ai/ Aq) quadratic > {Ai/ Ao)guartic- Note that 
the ratio of form factors also depend on the choice of the Fermi momentum of the B meson, 
as a smaller(larger) Fermi momentum would make the form factors more(less) sensitive to 
the tail of the wavefunction of M, as well as mixing effects in the wavefunction of the meson 
M. 

One can check that the predictions for the ratios of branching ratios are not very different 
from the case where we neglect radial mixings. Hence one concludes that the effects of radial 
mixing are in most cases small and negligible for practical applications. We observe in 
Table. |l| that there can be a large enhancement for R^. One can get a rough estimate of the 
branching ratio for Bg (pir^ from 



BR[Bs 07r°] _ 1 
BR[BO P+TT-] ^ 2 



0.04 



VubV:,{c2 + cl/N,) 

where we have neglected form factor and phase space differences between B^ (pn^ and 
B^ p+n-. Using the measured BR[B p+n-] ~ 28 x lO^^ [0 we get BR[Bs 07r°] ~ 
10^^. Hence the large enhancement for i?^ indicates that BR[Bs tp'n^] can be of O(10~^). 

Note that in the piu) system there are two resonances, p(1450)[cij(1420)] and p( 1700) [cj [1650], 
which can be identified a S-wave radial excitation(2S) and a D wave orbital excitation in 
the quark model. However recent studies of the decays of these resonances show that it 
is possible that these states are mixtures of qq and hybrid states RefJlOl- Hence the state 



p(1450)[a;(1420)] is interpreted as a 2S state with a small mixture of a hybrid state. We 
do not take into account such possible mixing with a hybrid state in our calculation and 
the meson masses for these excited states used in our calculation are the ones we predict in 
section 2. For the (j) system there is only state at 0(1680) which we interpret as a 2S state 
in the absence of mixing effects. 

We construct two ratios 

Rk = BR{B- K-7])/BR{B R-r]') (49) 

and 

Rk* = BR{B- K-*7]')/BR{B K-*7]) (50) 
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Let us assume the r] — rj' mixing used in Ref|]T7| 



1^) = 7f[^o-^o] 



V2 



[No + So] 



(51) 



where, as before, >= \uu + dd > /\/2 and S = \ss >. Now, from Eqn. |35| if we only 
include the term Pti then we find 



R 



K 



fKF+ + 



(52) 



where we have dropped the masses of the pseudoscalars in the final states and the form 
factor is defined through 



{P{Pf)\q7,il-l,)b\B{p,)) 



rrin — m. 



+ 



— m] 



-Qt^Fo 



Fp{pi+Pf)f, + Fpqf,. 



(53) 



In the above equation fx is the kaon decay constant and the decay constants and /^/ are 
defined by. 



(54) 



For the mixing in Eqn. ^ and assuming SU (3) flavor symmetry, we can write 

fu,d ^ Ik 



rs 
Jn 



-I 



K 



fr)' 



V2 

d _ fK_ 

2 
2 

V' 9 



(55) 
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One can then write 



R 



K 



1 l_ 

2 V2 



2 ^ V2 



0.03 (56) 



It was shown in Ref|T^, |T8[ that there is a parity selection rule in the decays B r](r]')K^*^ 
which fixes the relative phase between the penguin amplitudes representing the strange 
and nonstrange contributions to rj and rj' final states. In particular the parity selection 
rule predicts that the phases between the strange and nonstrange penguin amplitudes in 
B Kri{ri') and B K*ri{r]') are reversed. Hence neglecting form factor differences for 
B ^ P and B ^ V transitions, one obtains 



R 



K* 



1 3_ 

2 



2 ^ V2 



0.03 (57) 



We now calculate the ratios above with a nonstandard rj — rj' mixing. We will use the 
mixing for a quartic potential given in Eqn. This is because the ground states of this 
mixing has the largest overlap with the standard mixing in Eqn. |51|. Later on in this paper 
we will compare form factors calculated with the mixing in Eqn. ^ with those calculated 
with the mixing in Eqn. |51| to get an idea of the effects of radial mixing in the rj — rj' 
wavefunctions. A key ingredient in the parity selection rule that predicts the relative phases 
of the strange and nonstrange contributions to f] and r]' final states is approximate fiavor 
symmetry. In the fiavor symmetry limit the radial wave functions of vr, K, t] and f]' ( up to 
mixing factors in r]{ri') are all the same. One can then argue that since all are ground state 
wave functions with no nodes and a constant phase over the entire radial domain, fiavor 
symmetry breaking can change the radial shape and size of the wave function but will not 
reverse the phases. This is probably a reasonable assumption as long as we only have ground 
state wave functions. But when radially excited wave functions which have nodes are also 
considered, there is a phase reversal in the wave function. Note that with factorization, the 
B decay into the nonstrange part of the t] and rj' involves a point like form factor for the kaon 
and and a hadronic overlap integral for the ri{ri'). But the B decay into the strange part of 
the 7] and rj' involves a point like form factor for the r]{ri') and a hadronic overlap integral for 
the kaon. The penguin amplitude that involves a hadronic overlap integral for the T]{ri') can 
cause a phase ambiguity between the two penguin amplitudes because the wave function for 
the radial excitation changes phase and we are not able to exactly calculate fiavor breaking 
effects. 
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A simple model to simulate flavor symmetry breaking would be to define, for the rjlr]'), 
an effective wavefunction 



eff 



'N2 



(5^ 



where the Oj's are from Eqn. |1^ and ri and r2, representing flavor breaking effects, can have 
both signs. The effective wavefunction "^eff will then be used for form factor calculations. 
We will consider two choices of ri and r2. In the first case we set ri = r2 = 1 and in the 
second case we choose ri and r2 to be such that the contributions from the various radial 
excitations add constructively in the form factor calculations involving the rj and the rj'. 

The form factors in Eqn. 53 calculated using the quark model of Ref [^] are given by 



where 



A 

Jl 
J2 



4^. + ^^.] 
M 



AMiMfM+ 



Mi-q 

l3 



m 



d p<p*f{p + a)(t)i{p) 
J m TJif 



(59) 



a 



Mi±M 



f 



m^q 



,Ml-q' 



M 



■Ml-q^ 
rriimf 

rrii + TUf 



(60) 



and (pf and (pi represent the momentum space wave functions, (3 is the velocity of the mesons 
in the equal velocity frame and rrii f are the non spectator quark masses of the initial and 
final meson. For the calculation of the form factors in the ri{ri') system, clearly, the weak 
binding assumption, which would imply ~ M^/ ~ Mp, does not hold. One does not know 
how to incorporate corrections to the weak binding limit. One could choose, for instance, 
different quark masses, and there are various reasonable choices of quark masses that can be 
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made leading to very different predictions for tlie form factors. We do not wish to explore 
all these possibilities here, instead we will use the same procedure and the same values for 
the quark masses used in the calculation of the form factors in the vector system. We then 
find 



+nonstandard 
+standard 

+nonstandard 
+standard 



1.5,1.7 



0.5,2.1 



(61) 



where Fj^standard ^ud F^nonstandard ^rc the form factors calculated in the standard mixing in 
Eqn. and for the nonstandard mixing in Eqn. |l^. The contribution from the form factor 
F_ is negligible. The two numbers in the equation above correspond to the two choices 
for ri 2 mentioned above. We note that the form factor with nonstandard mixing does not 
change much for the rj' but changes significantly for the r] for the two choices of ri 2. 

Using the mixing in Eqn. ^ one obtains for the decay constants 



£u,d 

^nonstandard 

ru,d 

J ^standard 
rU,d 

In' 

I nonstandard 

ru,d 
In' 

I standard 
Vnonstandard 
^standard 



nonstandard 



fv's 



1 

1.1 
0.8 
1.2 



(62) 



We will choose the second entry in Eqn. BT] and for simplicity use 



+nonstandard ^ 

+standard 
ru,d,s ^ 

Vnonstandard 

ju,d,s ^ 

V nonstandard 



2.0 

ru,d,s 

Vstandar 



V standard 



+nonstandard 
+standard 



(63) 



The relations in Eqn. ^ are within the uncertainties in the calculation of the form factors 
in Eqn. |6l| and Eqn. ^ 

It is now easy to check that the predictions for Rk and Rk* in Eqn. ^ and Eqn. |57| 
remain essentially unchanged. However, as was shown by one of the authors of this paper. 
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one can use flavour topology characteristics of charmless B decays to derive additional sum 
rules connecting B decays to Kri{ri') and Kn final states. One of the interesting sum rule 
derived in [j^, |TB|, ^ is, neglecting phase space corrections, 



Note that this sum rule is true for any standard rj — rj' mixing. Recent experimental measure- 



ments [^, however, show the above sum rule to be invalid. In the factorization assumption. 



one can write, neglecting the tree contribution. 



R 





2 

+ 


JkFj^' + f^'Fx 




F^ofx 


2 



(65) 



With the mixing in Eqn. ^ we get R ^ 3. Note that one can check that for any standard 
1] — 7]' mixing one always gets R ^ 3. With the nonstandard rj — rj' mixing in Eqn. |T^, we 
get R ^ 6 which is now consistent with experiment. One would also get similar predictions 
with a K* in the final state. In particular we have 

As we have argued before the OZI suppressed terms may play an important role in the 
decays B rj{rj')K^*^ decays. In particular the OZI suppressed terms are important for 
decays with a rj' in the final state because they add constructively while for the rj in the final 
state the OZI suppressed terms tend to cancel among themselves. One can calculate the 
contribution to the amplitude of the OZI suppressed term with the mixing in Eqn. ^ as. 



1 1 

1 1 



x^, ^ -(! + —) -0.21 



X, ^ 0.037 (67) 

where we have dropped factors common to both Xn and x^'. We see that the OZI suppressed 
contribution for the rj' in the final state is indeed more important than for the rj in the final 
state. 

We now present the full amplitude for the decays B rj{rj')K{K*) including all the 
terms in the effective Hamiltonian in the factorization assumption and including the OZI 
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suppressed terms. To make definite predictions we choose the form factors, 



viv') 



+ standard^ 



as 



well as the form factors for B K(K*) transitions to be given by Ref[^. We will use the 
decay constants /^'''^ = ^ 0.8/, and f^, = -f^ ^ 1.3/, . 



Finally, one can write the amplitude for i? ^ K rj' as ||23[| 



where 



M 



72 



Vuia^nQK + a2Q^')- Yl V.. {(Tin + T^r2)QK + T^Q,'} 



i=u,c,t 



Ti 



Ti 



al + al- al + (2ag - < 



^ 2ms{mb - ms) 2^^ 2^^ 2^^° 



a\ + 2(al + al 



m 



K 



1 



m„ + msmb- m„ 



+ a 



10 



with ai = ci + a2 



C2 + jr, a, — c- H 



fU 

J„' 



r2 



(68) 



i-s 



j+1 ' N,^ ■ ^ - U ^ h 

Qk = iF(f (m2,)(m| - m^)/„ Q^. = iF^' {m]^){m% - m'^^,)fK, V^ = Vu,Vc,Vt and A^'^ is 
effective number of colors. 

In the above equations we have used the quark equations of motion to simplify certain 
matrix elements. The masses used in these equations of motion are the current quark masses 
given in section 3. The expression for the amplitude can also be used for B — > rjK by making 
the necessary changes. It is also straight forward to write down the amplitudes for B K*rj' 
and B — > K*ri decays. 



M 
where 
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Si 



SI 



Vu{air;,A + a2mK*gK'B)- V,{{S\f:,> + S\f:^,)A + SlruK^gK^B] 



i=u,c,t 



1 



1 



24 - 24 - + 2^9 
03 + 04 - - {2al - al) 

~|~ Cl]^Q 



mt 



2ms{mb + rus 



+ o«7 



no 



(69) 



with A = 2mK*Aoe* ■ ps, B = 2e* ■ psFi {jn\,) and we will use qk* =221 MeV where 
gx* is the vector meson decay constant. A similar expression for B —>■ r]K* can also be 
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obtained. To identify the various terms in Eqn. ^ and Eqn. ^ , let us for simplicity drop 
the electroweak penguin terms represented by a\ — a\Q. Note that the term proportional to 
is formally of 0{l/mb) and so strictly in the rrib —>■ oo limit this term vanishes. However 
for realistic quark masses this term is not negligible and is chirally enhanced because of 
the strange quark mass, rris which is given in section 3. The OZI suppressed terms are 
represented by the terms a\ and a^. One can check that if we drop the OZI suppressed 
terms as well as the chirally enhanced terms then, neglecting form factor and phase space 
differences, we would recover the predictions of Eqn. |5^ and Eqn. |5^ from Eqn. ^ and 
Eqn. If we include the OZI suppressed terms in the calculation without the chirally 
enhanced contribution we obtain BR[B rj'K*] = 1.1 x 10"^ while if we ignore the OZI 
suppressed terms then we obtain BR[B — » r]'K*] = 0.34 x 10^^. Hence we see that the 
presence of the OZI suppressed term can alter significantly the BR[B rj'K*]. There is 
a much smaller effect of the OZI suppressed terms in BR[B v'K]- This, as already 
mentioned, is due to the fact that the OZI allowed terms in BR[B —>■ rj'K] tend to add 
constructively while they add destructively in BR[B rj'K*] and so the OZI suppressed 
effects are felt more strongly in B ^ rj'K*. Finally we present in the Table. ^ the results of 



Table 2: Branching ratios(BR) for B — > ri{ri')K{K*) decays 



Process 


Experimental BR 24 


Theory BR 


B- K-T]' 


{80+1° ± 7) X 10-6 


93 X 10-6 


B- K-7] 


< 6.9 X 10-6 


1.04 X 10-6 


B- K-*ri 


< 35 X 10-6 


3.6 X 10-6 


B- K-*r] 


(26.4^1^ ± 3.3) X 10-6 


10 X 10-6 



our calculation including all the terms in Eqn. |68| and Eqn. We see from Table. |^ that 
our calculations are in reasonable agreement with experiment. In particular we note that the 
chirally enhanced contributions tend to further increase the branching ratio BR[B rj'K*]. 
From the table above we can calculate 

Rk = 0.01 

Rj^, = 0.36 (70) 

If experiments find that the ratio Rk* is indeed much smaller than predicted here then this 
would indicate the presence of large non factorizable corrections that effectively cancel the 
OZI contributions as well as the chirally enhanced corrections. 
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6 Summary 



We have considered the weak decays of a B meson to final states that are mixtures of 
S-wave radially excited components. We calculated nonleptonic decays of the type B 
P'-k/B pn, B u'-k/B ujTi and B —>■ (P'tt/B —>■ (pn where p', uj' and (j)' are higher 
p, UJ and (j) resonances. We found that the transitions to the excited states can be com- 
parable or enhanced relative to transitions to the ground state. It would, therefore, be 
extremely interesting to test these predictions. We also studied the effect of radial mixing 
in the vector and the pseudoscalar systems generated from hyperfine interaction and the 
annihilation term. We found the effects of radial mixing to be small and generally negligible 
for all practical purposes in the vector system. However, in the rj — rj' system the effects 
of radial mixing are appreciable and seriously affect decay branching ratios. In particular 
we found that the experimental violation of the sum rule for B —>■ Kr]' in Eqn. |6^ can 
be explained by radial mixing without need for the OZI suppressed transitions. We also 
pointed out that the place to look for an OZI suppressed contribution is in 5 — »• K*!]' decays 
where the the OZI suppressed transitions become important as OZI allowed term is small. 
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